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We introduce the notion of a local shadow for a black hole and determine its shape for the particu¬ 
lar case of a distorted Schwarzschild black hole. Considering the lowest-order even and odd multiple 
moments, we compute the relation between the deformations of the shadow of a Schwarzschild black 
hole and the distortion multiple moments. For the range of values of multiple moments that we 
consider, the horizon is deformed much less than its corresponding shadow, suggesting the horizon 
is more ‘rigid’. Quite unexpectedly we hnd that a prolate distortion of the horizon gives rise to an 
oblate distortion of the shadow, and vice-versa. 
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I. INTRODUCTION 

The mysterious radio source, called Sagittarius A*, 
at the ceutre of Milky Way, with a mass of arouud 
4 X 10^ Mq is cousidered to be a supermassive black 
hole. To a distaut observer, the black hole casts a rela¬ 
tively large “shadow” with au appareut diameter of ^ 10 
gravitatioual radi|^ due to beudiug of light by the black 
hole. The predicted size of the shadow of Sagittarius A*, 
30/iarcsec), will be probed iu the uext few decades, 
usiug the Eveut Horizou Telescope via Very Loug Base- 
liue luterferometry [1]. 

The origiual iuvestigatiou of the shadow of a 
Schwarzschild black hole [2] has beeu exteuded to a uum- 
ber of other cases, iucludiug Kerr black holes [3], (see 
also mi), Kerr-Newmau black holes [7], Schwarzschild- 
de Sitter black holes [8], Tomimatsu-Sato space-times [9], 
black holes iu exteuded Cheru-Simous modified gravity 
m, rotatiug braueworld black holes mi ng, Kaluza- 
Kleiu rotatiug dilatou black holes [13], for the Kerr- 
NUT black holes [14], Kerr-Newmau-NUT black holes 
with a cosmological coustaut mi, aud multi-black holes 
[T6] . Shadows for black holes iu modified gravity m, 
aud the five-dimeusioual rotatiug Myers-Perry black hole 
m have also beeu iuvestigated. The effect of scatteriug 
ou the visibility of the shadow was uumerically demou- 
strated by Falcke, Melia aud Agol m- The visual ap- 
pearauce of au accretiou disk was studied with the help 
of various ray-traciug programs by several authors [20]- 
[23] . The overview of observatious aud the simulatious of 
pheuomeua for the black hole iu the ceutre of our galaxy 
was studied iu m- 

Most of these iuvestigatious of black hole shadows are 
for idealized cases, dealiug with exact solutious iu the ab- 
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^ 'i’he diameter of the “shadow” will be slightly different in the 
case of a rotating black hole. 


seuce of matter. It is clearly desirable to seek alteruatives 
to this case, particularly those for which exterual sources 
of matter iuflueuce the black hole. ludeed, to uuderstaud 
the properties of black holes as astrophysical objects pre¬ 
dicted by the theory of geueral relativity, it is esseutial to 
study their iuteractiou with exterual matter sources. A 
particularly relevaut example from astrophysics is that of 
a black hole iu a biuary system, where tidal effects from 
its partuer distort the black hole. Although the best way 
of studyiug such dyuamical systems is uumerical aual- 
ysis, obtaiuiug the most geueral dyuamical exact black 
hole solutiou iuteractiug with a geueral matter coufigura- 
tiou is out of reach. However cousiderable iusight cau be 
gaiued from studyiug exact solutious describiug a black 
hole distorted by exterual matter. Distorted black holes 
have beeu of cousiderable iuterest for quite some time 
[25H34] . Such solutious approximate those of dyuamical 
black holes that relax ou a time scale much shorter thau 
that of the exterual matter. 

Here we iuvestigate how the shadow of a distorted 
black hole gets deformed. The space-time we cousider 
uot asymptotically fiat because the sources of distortiou 
are exterior to the black hole regiou we study: the ob¬ 
server of the shadow is uot located at asymptotic iufiuity, 
but rather at a fiuite distauce from the black hole, iu a 
vacuum regiou that is iuterior to the exterual sources. 
We are iuterested iu situatious where the distortiug po- 
teutial does uot domiuate over the poteutial of the black 
hole. For this reasou we restrict ourselves to the study 
of small distortious, aud for simplicity we cousider ouly 
the simplest such types. We defiue the coucept of a ‘lo¬ 
cal shadow’, a shadow seeu by au observer uot located 
at asymptotic iufiuit30 Oue of our key results is that of 
relatiug the deformatious of the local shadow of the black 
hole to its distortiou parameters such as the quadruple 
aud octupole multiple momeut of the distortiou fields. 
We also fiud eyebrow structures for a siugle distorted 


^ The shadow for an observer not at infinity has been considered 
in HU- However, our definition of local shadow is different. 
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black hole, similar to that which occurs for two merging 
black holes [35]. 

In Section 2, we first review the solution representing 
a distorted axisymmetric Schwarzschild black hole. In 
Section 3, we construct the null geodesic equations for 
this space-time. In Section 4, we construct a map which 
would enable us to define the local shadow. In Section 5, 
we present and discuss our results, and finish with Section 
6 where, we sum up. In what follows employ units where 
G = c= 1. 


II. METRIC 

We consider a Schwarzschild black hole in an external 
gravitational field. The solution is static and axisymmet¬ 
ric and its metric can be presented in the following form 

m, 


dS^ = - (l - 

-h dO^ sin^ Odc /)^), (1) 

Here, U = U{p^ 0) and V = H(p, 0) are functions of p and 
0. For U = V = 0, this metric represents a Schwarzschild 
black hole. The null rays in this space-time have four 
conserved quantities, energy f, angular momentum L, 
the azimuthal angular momentum and = 0, 

where = dx^ jdr is the photon 4-velocity, with r an 
afhne parameter along a null geodesic. It is convenient 
to write the metric Eq. 0 in dimensionless form, 

where 

dS‘^=rldS‘^, T=-, p = rgr, (3) 

Vg 

with Tg = 2M. The explicit form of U and E as a func¬ 
tions of r and 0 are 

U=J2cnR^Pn 

n>0 

n —1 

V = [cos^ — (2r — 1) 

n>l z>o 

-(-l)"-*((2r-l)+cos0)] RXi 

+ E '^^R^^'^iPnPk-Pn-lPk-l) (4) 

n ^ k 


where, Pn are the Legendre polynomials of the first kincQ 

P„ = P„((2r-l)^) 

R = Y^(27^^-T)^'^-shi^. (5) 

Moreover, the metric can be written in the form 

-1 

ds2 _ g-2no 


_(1 _ + (1 - b e-^+^^dr^ 


+ e ^^r^{e^'^dd‘^ + sin^ 0d(p^)] , 


(6) 


where 


U = U-uo, = t = Te^^° (7) 


n>0 


In order to avoid the conical singularities in the space- 
time we must have 

C2n+1 = 0. (8) 

n>0 


In what follows, we construct the null geodesic equations 
of the following metric 

ds^ = -(1 - -)e^^dt^ + (1 - l)”b-2^^+2V'd^2 
r r 

+e-^^r^ie^^de^ + sin^ Odcp^), (9) 

Note, that metric Eq. 0 is related to Eq. 0 by conformal 
transformations and redefinitions of the time coordinate. 
The null geodesic equations are conformally invariant. 


III. NULL GEODESICS 

To construct the shape of the shadow of the space-time 
Eq.0, we need to study the null geodesics for this space- 
time. As with Schwarzschild space-time, the space-time 
0 possesses two commuting Killing vectors 

= (10) 

corresponding to the symmetries of time translation and 
rotation around the ^ axis. The corresponding quantities 
conserved along geodesics of the space-time are 

Lz = sin^ 0 if) (11) 


^ Legendre polynomials are solutions to Legendre’s differential 
equation, 


dx 


(1 - X^)—Pn{x) 
dx 


+ n{n + l)Pn{x) = 0 . 


They can be written by using Rodrigues’ formula, 


Pn{x) = 


d^ 


2'^n\ dx'^ 


[(x2 - 1)-] . 
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respectively, the photon energy E and azimuthal angu¬ 
lar momentum Lz- Here = dx^jdr is the photon 
4-velocity, normalized as = 0, with r an affine pa¬ 
rameter along a null geodesic. The geodesic equation is 


d?x^ dx^ dx^ 

-=-h - 

dr dr 


( 12 ) 


are the Christoffel symbols defined in terms of the 
metric tensor. From Eq. we see that d^/dr has 
monotonic behaviour, never changing sign along a par¬ 
ticular trajectory. 

Therefore, we introduce 


dr dr dcj) , 

dr dcj) dr 

dO dO dcj) 

dr dcj) dr 


(13) 


Here, the overdot denotes the derivative with respect to 
the affine parameter, and prime defines the derivate with 
respect to (j). The geodesic equations take the following 
form 


^/2 




hr r 3 — 4r 

+ 2 ^ H— ? -TT 

h j r[r — 1) 


nf2 

+ ^(r-l) 


hr 

+ 2 

h 


n/ f h o ^ « 

Or - - - 7 —h 2 cot 6 

f h 


9” = 


sin^ 9 

2llhp{r - 1 ) 

T sin^ 9{r - 1) 

r'^h^o 9'^ 
2 r(r - l)h ^ T 


l + r(r-l) 


f,r 


rf,r 

2 / 


(14) 


h,e , A A-Q 

--7 - 2 ^ +4cot6» 

h f 


h^r f^r 

~h !. 


9'r' 

„;„2, 


sin^ 9f^o sin^ 9 


2llhnr - 1) 2hf 


E — 2 cot 9 

2 / 


(15) 


where, / = h = , r" = d?r/d4>^, 9” = 

cP9/d(j)‘^, and we have defined Iz = LzjE. We also have 
the following constraint 


other n. The other case is the quadrupole moment where 
Co 7 ^ 0 and C 2 7 ^ 0 and c„ = 0 for any other case. We 
shall consider the dipole and quadrupole cases separately. 
Recall that since we made the conformal transformation 
the quantity Cq does not contribute to our result. 

Here, we have assumed that in analogy to the Newto¬ 
nian multiple moments we have cq > ci > C 2 > ... > 
Note that if we calculate the multiple moments for a 
static axisymmetric Newtonian gravitational potential of 
two point-like masses Mi and M 2 and a ring located 
on the equator of the axis joining the two masses, we 
find that Cq > Ci > C 2 > ... > c^, and the ring con¬ 
tributes only to the even multiple moments. In this sense 
a quadrupole multiple moment C 2 is analogous to that of 
a ring around the black hole. For the quadrupole distor¬ 
tion we have 


V = — 2 c 2 sin^ 0 [2r — 1 18 c 2 r^ cos^ 0 — 36c2r^ cos^ 0 

+ 22c2r^ cos^ 0 — 4c2r cos^ 0 — 2c2r^ 

+ 4c2r^ - 2c2r^] , (17) 

U = —C 2 ^{2r — 1)^(3 cos^ 0 — 1) ^ sin^ 0. (18) 

For the dipole distortions (odd multiple moments) we 
have 


V = A cos 0 ^ Aq cos^ 0 ^ Ai cos^ 0 ^ A 2 
hi = —Cl cos6> [r [cos^ 6>(20r^ — 30r -j- 12) 

- (12r^ - 18r + 8) ] + sin^ (9] (19) 

where the functions A, Hq, Hi, H2 of (r,d) are given by 

A = 2ci sin^ d [6r(r — 1) + 1] , 

Ho = -2ciVsin^ (9 [-168r^ + b04r^ - 570r^ 

+ 300r^ - 72r + 6] , 

Hi = —2ci^r sin^ 0 [300r^ — 900r^ + 1008r^ 

- 516r^ + 117r-9] , 

A 2 = - 2 ciVsin^ 0 [I 2 r^ - 36r'^ + 42r^ 

- 24r^ + 7r - 1] . (20) 


IV. MAP 


= 0 , 

sin^ OF + fhr'^ + r^FfhO'^ = 0 

f H 

(16) 

where F = 1 — 1 /r. 

In order to study the effect of distortions on the black 
hole, we consider the simplest cases. The simplest type 
of distortion is due to a monopole cq 7 ^ 0 and c^ = 0 
for n > 0. The next less trivial distortion is due to a 
dipole Co 7 ^ 0 and ci 7 ^ 0. In this case from condition 
(§, we have ci = — C 3 , and we consider c^ = 0 for any 


Our goal is to see how an observer located at (ro,^o) 
sees the shadow of a distorted black hole. Consider rays 
emitted from any point on a sphere of radius Ce around 
the black hole that may reach the observer located at 
i'^o^Oo) or else be absorbed by the black hole as shown 
in figure This observer is located in the interior re¬ 
gion, which means that the sources responsible for the 
distortion of the black hole are all located at r > > Ce- 

To analyze the problem we consider rays emitted at 
(ro, 6 >o) and trace them backward to their source. In or¬ 
der to characterize properties of rays emitted at (ro^Oo), 
which are directly connected with observation, we pro- 
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ceed as follows. First we choose the following orthonor¬ 
mal tetrad 




( 21 ) 


at the point of the observation, where ... 


stands for 


the limit where r = Vq and 0 = Oq. These vectors are or¬ 
thonormal. In what follows, we use (t, r, 6>, 0) coordinates 
and the indices are (0,1, 2, 3), respectively. The vector ei 
is in radial direction to the direction of the black hole. 
The tangent vector to a null ray is 


u 


M _ 


dx^ 

dr 


( 22 ) 


We project this vector into the observer’s orthonormal 
frame. Thus, the tangent vector at the observation point 
(^0 7 ^o) can be written as 

«'^ = /3(-e^ + C'e5‘ + CM+CM), (23) 


where and are displacement angles. Replacing 


11 


Eq.(^ in Eq.(|23[) and comparing to Eq.(|22[) and using 
Eq.( 


, we find the value of the scalar coefficient P 


P = - 


E 




(24) 


Consider a photon that gets emitted from a point lo¬ 
cated on a sphere of radius rg. Eollowing the motion of 
this null ray forward, it will either get captured by the 
black hole or reach the eye of the observer located at 
^o) depending on the shape of the black hole. To de¬ 
termine the shadow of the black hole, we instead trace 
the trajectory of these photons backward (see figure [^. 
The quantities give the angle of the photon as 

it reaches the observer located at point (rg, ^o). Suppose 
the observer is looking in the direction of the centre of the 
black hole. Tracing straight backward null rays leaving 
the observer at angles as if the space were fiat, it 

appears to him/her that the photon has reached his/her 
eye from the point and on the plane 

of the black hole. If this null ray is absorbed by the black 
hole, we consider to be a black point on the plane 

of the black hole - by definition it is a member of the 
“local shadow” of the black hole, so named because the 
observer is not located at asymptotic infinity. Replacing 
Eq.([2T]) in Eq.([2^ and comparing to Eq.(22), we find the 
displacement angles , and and using Eq. (24). More¬ 


over, by multiplying by and using Eq.(13) we find 


f =±-T 

f = - 


sin^ 0 
I 


fhiHi-lw 


sm 6 \ r 


(25) 


f Z 


► y 



FIG. 1. The photon emitted from a point on the sphere of 
radius Ve reaches an observer located at (ro,^o)- Ge and 
are unit vectors of the observer’s orthonormal tetrad directed 
along the coordinate lines of 0 and 0. The red dashed line 
represents an example of a null trajectory that reaches the 
observer. 


On the other hand, if this null ray reaches a radius rg, 
(after propagating in the space-time) it is not a member 
of the local shadow of the black hole. 

Before we close this section of the map we perform a 
“test”. We want to see if by taking the limit r ^ oc and 
setting the distortions to zero (/ = h = 1) we will get 
the map [3], 

- • (26) 

r^oo 

At infinity and without distortions the system is com¬ 
pletely integrable and the equation of motion for the 0 is 
given in the following form. 



0 = 


y/Lo^ sin^ e - L/ 
sin 0 


Therefore Eq.(26) takes the form, 


f = 


h 


sin^ 6 


For we have, 


f = -w 


sin6> V r 


sinO 


(27) 


(28) 


(29) 


We see that equations ( 28|29 ) match with the equa¬ 
tions for the map of undistorted black holes [3]. In the 
Schwarzschild undistorted space-time, Eq. 0 for U = 
E = 0, there exists a second rank Killing tensor 

= 0) which generates another conserved quan¬ 
tity, the squared total angular momentum of a photon 


P = 




(30) 
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For the critical impact parameter ^critical ~ j = 
27M^, the light enters a circular knife-edge orbit of radius 
p = 3M, and may orbit the black hole for part of an orbit 
or for many turns before it escapes or plunges. These rays 
barely escape the Schwarzschild black hole and form the 
rim of the black hole shadow. For the rescaled metric 
Eq.([^ with U = V = These rays have L‘^ jE‘^ = 27/4 
where, L = Ljrg and £ = VgE. For an observer located 
at infinity this implies that 

(?)' + (?)’ = f. (31) 

Note that, usually the radius of the shadow of a 
Schwarzschild black hole as seen by an observer at infin¬ 
ity is 3\/3M. However, in our case, since, p = 2Mr, (see 
Eq.([^), the radius of the shadow of a Schwarzschild black 
hole as seen by an observer at infinity is \/^/4 = 3\/3/2. 
Erom equation ( [^ , when / = h = 1, for an observer lo¬ 
cated at To, we have 


(^-2)2 ^ (^-3)2 


Iz’^jro - 1 ) 

sin^ 9ro 


(6>'2 + sin^ e) 


27{ro - 1 ) 



upon setting = 27/4 and using = 

sin sin^ 0 — Iz^ from Schwarzschild space-time. 


V. RESULTS AND DISCUSSION 


Here we study the effect of distortions on the local 
shadow of an undistorted black hole. Imagine for exam¬ 
ple that there is a ring around the black hole with radius 
R, and the observer is at a radius r^ < R. The distortions 
are characterized in terms of the multiple moments c^. 
Eor simplicity we consider the lowest order multiple mo¬ 
ments (considering odd and even moments separately). 
In choosing the appropriate value of multiple moment, 
we suppose its value to be small enough that the ratio 
/ = gttd/gtt of the tt-component of the distorted black 
hole metric to its Schwarzschild counterpart is not very 
large, i.e., / < 10. This ratio is illustrated in Eigurej^ 


We have also calculated the / and h functions in figures 
We need to consider smaller values of ci compared to 
C 2 . Eor example, ci = 1/150 is a very strong distortion, 
such that gttdigtt > 10000. We need to keep in mind that 
we want to consider a black hole distorted by external 
sources that are not much stronger gravitationally than 
the black hole itself. We consider —1/150 < C 2 < 1/150; 
for this range of values of C 2 we have, gttdigtt < 2.6. We 
assume, —1/800 < ci < 1/800; we find gttd/gtt < 6 for 
this range of values of ci = — C 3 . 

In the following, we consider only the case where 
Te = Vq. In our numerical computations, for every photon 
trajectory with initial values and we compute Iz, 
O'. Erom the constraint equation Eq.(16) we can derive 


r' as following 

sin^ 0 


fHlh 




fh 


(33) 


Erom Eq.(33), we see that and cannot take every 


value. Computing this quantity for any choice of {vo^Oq) 
yields 


r\ro,0o) = 


sin^ 0 


pllh 


(ro,0o) 


(34) 


and for and outside a circle of radius r^, the ex¬ 
pression under the root becomes negative. Thus relation 
Eq.(34) gives us the range of values of and that are 


valid. Numerically, we divide the interior region of the 
circle of Eq.(34) to small pixels of side 0.01; for each of 


these values we compute the corresponding photon tra¬ 
jectory. We see whether this trajectory ends inside the 
black hole or escape the black hole, namely reaches ra¬ 
dius To- Eor each trajectory, we calculate the maximum 
deviation from the constraint Eq.(16). Erom this we find 


the maximum deviation from the constraint for all the 
trajectories, which we find to be 4 x 10“^. 

We present the effect of distortions for even multiple 
moments in figures [5][6l We see that for an observer 
located in the equatorial plane, Oq = 7r/2 the shadow 
looks like an ellipse. Eor positive and negative even 
multiple moments C 2 , we respectively obtain prolate and 
oblate shadow shapes. As the multiple moment increases 
through positive values, the shadow becomes more pro¬ 
late, whereas it becomes more oblate for increasingly neg¬ 
ative values. Eor an observer not located on the equato¬ 
rial plane, the local shadow of the black hole moves up¬ 
ward or downward in the impact plane of observer. As 
the magnitude of the negative multiple moment increases, 
the shadow moves down (when is the horizontal axes, 
illustrated in figure . As the value of positive multi¬ 
ple moment increases the shadow moves up, as shown in 
figure . 

Since the shadow looks like an ellipse for an observer 
located at = 7 r/ 2 , we fit the shape of the shadow to an 
ellipse, taking a and b to be its semi-major/semi-minor 
axes, the origin at the centre, and the angular coordinate 
(/) to be measured from the axis To find a and b we 
minimize the normalized square error 


J = 




Pi? 


EZoip~? 


with 


(42) 


pi — 


- =[(f»)“+(f 

ab 

cos^ 4)1 + a? sin^ c, 


(j)i = arctan(^), 


(43) 


where is given by our numerical result, and pi is 

the equation of the ellipse for every (pi. In figure|^we plot 
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FIG. 2. Comparison of gud of the distorted black hole FIG. 3. / function for distorted black hole. Top: ci = 1/800; 
(black), to gtt of the Schwarzschild black hole (grey). Top: Bottom C 2 = 1/150 

Cl = 1/800; bottom for C 2 = 1/150. 


the parameters a and 6 as a function of the quadrupole 
moment C 2 for an observer located on the equator {Oo = 
7r/2) and at a fixed distance Vq from the centre of the 
black hole. To find this relationship, we assumed that 
a and b are linear functions of C 2 and that the intercept 
is given by the value (27(ro — l)/4ro)^/^ (see equation 
)). Obtaining the slope, we see that the assumption 
of linearity between a or 6 and C 2 is a good one, since it 
gives a fitting error comparable or smaller than the size 
of our pixel, (recall that this is 0.01). For the observer 
on the equatorial plane and distance To = 4, we find 


(32 


a = -99.68C2 + 2.25, b = 22.58c2 + 2.25, (45) 

where the respective errors in the fits for (a, b) are (4.3 x 
10 - 1.1 X 10 ^). As stated, these errors are smaller than 

the size 0.01 of our pixel, and so this is the actual error. 
For C 2 = 0 the shadow is a circle of radius 2.25. For an 


observer on the equatorial plane at To = 5, the values are 


a = -161.96C2 + 2.32, b = 34.11c2 + 2.32, (46) 

where the fitting errors for (a, b) are now (3.5 x 10“^, 1.2 x 
10“^) respectively. For C 2 = 0 the shadow is a circle of 
radius 2.32. As the observer moves away from the black 
hole the radius of the circle increases. For the undis¬ 
torted Schwarzschild geometry (c 2 = 0 ), and the shadow 
is a circle of radius 3V3/2 ^ 2.598 (see Eq.(31)) for an 
observer located at infinity. Table |T] indicates that as the 
observer is moves away from the black hole the magni¬ 
tudes of the slopes of the parameters a and b increase. 
This is also valid for the intercepts. 

Though it seems that the parameters a and b vanish 
for large values of C 2 , we emphasize that these slopes 
are valid only for small values —1/150 < C 2 < 1/150 of 
the quadrupole moment. Had we included larger values. 
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TABLE L Obtained values and errors of the parameters a and b for an observer on the equatorial plane at Vq. 


To 

a 

error in a htting 

b 

error in b htting 

3 

a = -41.85C2+ 2.12 

1.34 X 10"^ 

5 = 8.89c2 + 2.12 

7.25 X 10"^ 

4 

a = — 99.68c2 + 2.25 

4.3 X 10“^ 

b = 22.58C2 + 2.25 

1.1 X 10“® 

5 

a = — 161.96c2 + 2.32 

3.5 X 10“^ 

5 = 34.11C2 + 2.32 

1.2 X 10“^ 

6 

a = — 251.63c2 + 2.37 

2.4 X 10“^ 

b = 44.27C2 + 2.37 

3.3 X 10“® 



-4 -3 




(37) 

FIG. 4. h function for distorted black hole. Top: ci = 1/800; 
Bottom: C 2 = 1/150. 


say —1/5 < C 2 < 1/5, then the slopes in Table would 
have been totally different, since distortions drastically 
dominate over the black hole potential for larger values 
of the multiple moments. 

The comparison of the shadow of a distorted black hole 
to the embedding of its horizon is interesting. The hori¬ 
zon surface, t = const, r = 1, of the distorted black hole 
Eq. is given by 

dSjj = + sin^ 0d<f) (47) 


- 4 J (38) 

FIG. 5. On the top we can see the shadow of a black hole for 
an observer at = 7r/4 and radius ro = 5; the values of the 
multiple moments are C 2 = (with red, (3)) C 2 = 

(with black, (2)) and C 2 = 0, (which is the undistorted case 
with blue, (1)). On the bottom we can see the shadow of a 
black hole for an observer at Oo = tv/2 and radius ro = 5; 
the values of the multiple moments are C 2 = — , (with red, 

(3)) C 2 = (with black, (2)) and C 2 = 0, (which is the 

undistorted case with blue, (1)). 


where, Uh = = 1,0) and Vh = V{r = 1,0). Con¬ 

sider isometrically embedding this 2-dimensional axisym- 
metric metric into the 3-dimensional space 

= edz^ + dp^ + (48) 

where {Z, p, (j)) are cylindrical coordinates. Setting e = 1 
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-3 


(39) 




(40) 


FIG. 6. On the top we can see the shadow of a black hole 
for an observer at = Tr/d and radius ro = 5; the values of 
the multiple moments are C 2 = (with red, (3)) C 2 = 
(with black, (2)) and C 2 = 0, (which is the undistorted case 
with blue, (1)). On the bottom we can see the shadow of a 
black hole for an observer at Oo = 7r/2 and radius ro = 5; 
the values of the multiple moments are C 2 = p|q, (with red, 
(3)) C 2 = (with black, (2)) and C 2 = 0, (which is the 
undistorted case with blue, (1)). 


FIG. 7. On the top we can see the shadow of a black hole 
for an observer at = 7r/4 and radius ro = 5; the values of 
the multiple moments are ci = (with red, (3)) ci = 

(with black, (2)) and ci = 0, (which is the undistorted case 
with blue, (1)). On the bottom we can see the shadow of a 
black hole for an observer at 9o = tv/2 and radius ro = 5; 
the values of the multiple moments are ci = (with red, 

(2)) Cl = (with black, (3)) and ci = 0, (which is the 
undistorted case with blue, (1)). 


corresponds to Euclidean space, whereas e = — 1 corre¬ 
sponds to pseudo-Euclidean space. We have Z = Z{0), 
and p = p{0). Matching the metrics ( [47| ) and ( [4^ , we 
derive the following embedding map 

p{0) = sin 6 >, 

- pV)] ^ (49) 

J0 

Eor even multiple moments the horizon is deformed to an 
ellipse. Using the same fitting method as before, we find 

a = 1.07 c2 + 1, h= -2.11C2 + 1 (50) 

where we have fixed the intercept of the linear fit to be 1 
and found the slope; the errors in the above fittings are 
of order 10 “^. 


Quite counterintuitively, for positive C 2 the horizon of 
the black hole is oblate whilst its shadow is prolate. Eor 
the negative C 2 the opposite occurs: a prolate horizon 
has an oblate shadow. As noted above, a positive C 2 
corresponds to a ring on the equatorial plane in Newto¬ 
nian gravity. We can see from the bottom plot of fig¬ 
ure that for C 2 > 0, gud of the distorted black hole 
around the equator is less than that of an undistorted 
Schwarzschild black hole, whereas on axis gud is greater. 
This implies that the gravitational field on the equator 
for a distorted black hole is greater than the undistorted 
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FIG. 8. On the top we can see the shadow of a black hole for 
an observer at and radius ro = 5; the values of the 

multiple moments are ci = (with red, (2)) a = 

(with black, (3)) and ci = 0, (which is the undistorted case 
with blue, (1)). On the bottom we can see the shadow of a 
black hole for an observer at Oo = 7r/2 and radius ro = 5; 
the values of the multiple moments are ci = (with red, 

(3)) Cl = (with black, (2)) and ci = 0, (which is the 

undistorted case with blue, (1)). 


case and weaker on axiQ On the equator there exists the 
effect of the gravitational field of both the ring and the 
black hole. For a quadrupole moment C 2 the equations 
of motion imply that on the equatorial plane ^ = 0 when 
^ = 0 (note that and h^o vanish on the equator in this 
case). This implies the existence of planar motion. Since 
the equations of motion are integrable on the equator, 
we can construct an effective potential V for the planar 


Note that for C 2 > 0 the function 


U = - 


2 


2 


i2r-lf , 


on the equator is always negative. Therefore / = gttdidtt < 1- 


FIG. 9. The parameters of the ellipse a (red, (2)) and b 
(blue, (1)) vs. the multiple moment C 2 . Top: observer at 
Vo = 4; Bottom: observer at Vo = 5. 


motion of null rays m 


where 


1(^)2 = l _ y 2 

iVdt’ n 


df=[/(l-i)]V2, 


(51) 


(52) 


The maximum of the effective potential is equal to the 
value of 1//^ or 1/^criticaZ knife-edge orbit. For 


C 2 = 1/150 we have = 0.1389. Using Eq.(|25[), 

for an observer located on the equator and = 5 this 
corresponds to = 1.3894, which shows that the shadow 
of the distorted black hole on the equator has a smaller 
diameter than that of an undistorted Schwarzschild black 
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FIG. 10. ^3 as a function of the position of the observer in 

the equatorial plane.The blue, (1), colour corresponds to the 
undistorted case whereas the black, (2), colour corresponds 
to the distorted one. The value of C 2 = 1/150. 


hole, for which = 2.32 confirming the result of a pro¬ 
late shadow for C 2 > 0. The knife-edge orbit is located at 
radius r = 1.50, (or, p = 3M), for the undistorted case. 
For C 2 = 1/150 this orbit is at r = 1.46. Evidently the at¬ 
tractive influence of the ring allows potentially escaping 
null rays to more closely approach the black hole. 


The angular momentum of the photons for the knife 
edge orbit is larger than the undistorted case. If instead 
of Eq.(25), which is valid for a local observer, we had 
the map Eq.([^ (which is for an observer at infinity) 


this would immediately imply that ^3 would be larger for 
a distorted black hole. However the function / plays an 
important role in Eq.(25); on the equator it is always less 
than one and decreases as r increases. As a consequence 
^3 decreases with distance, implying the horizontal axis 
of the ellipse decreases in size as the observer gets further 
and further away from the black hole. Eor an observer 
very close to the black hole ^3 has almost the same value 
of that of a Schwarzschild black hole. This is illustrated 
in figure pT| 


While our metric represents a vacuum space-time, 
there must be (more distant) matter sources causing dis¬ 
tortion of the black hole. Provided these sources are lo¬ 
cated in some finite-sized region we can extend the so¬ 
lution beyond this (non-vacuum) region to a yet more 
distant vacuum region. In this exterior region we must 
replace the expansions of the functions U and V with ex¬ 
pansions in terms of exterior multiple moments. In this 
exterior vacuum region the function / increases until be¬ 
comes unity at infinity. Assuming the angular momen¬ 
tum of the photon is conserved as it (traced backward) 
passes through the non-vacuum and exterior vacuum re¬ 
gions, an observer in the exterior vacuum region can use 
the critical value of the angular momentum of the pho¬ 
ton to deduce ^ 3 . Eor an observer located at infinity ^3 is 


larger than that of the undistorted Schwarzschild black 
hole. 

In figure pH] we illustrate a schematic representation of 
an oblate horizon and a prolate shadow. All the observers 
in the (x, y) plane (at any 0 angle) see the same shadow 
shape. Eor illustration purposes the projection of the 
shadow was made on the other side of the black hole 
rather than on the plane of the black hole. One can say 



FIG. 11. Schematic representation of an oblate horizon and 
a prolate shadow. 


that the black hole horizon is more “rigid” than that of 
the black hole shadow, i.e. a larger value of the multipole 
moment C 2 is required in order to deform the horizon as 
much as the black hole shadow. 

Note that for very large values of multiple moments 
(which we do not consider), the horizon cannot be em¬ 
bedded in Euclidian space. Eor a 2-dimensional axisym- 
metric metric if the Gauss curvature is negative at the 
fixed points of the rotation group, it is impossible to iso- 
metrically embed a region containing such a fixed point 
in Euclidean space E^. Such surfaces can be globally em¬ 
bedded in E^ [36]. In the case of Schwarzschild black 
hole, with quadrupole distortion the Gauss curvature for 
C 2 > 1/12 becomes negative at both of the poles, ^ = 0 , 
and 0 = IT. Similarly, for the octupole distortion, the 
Gauss curvature becomes negative at one of the poles 
for |c 3 | > 1/20. Eor multipole moments greater than or 
equal to these values, the horizon surface of the distorted 
black hole cannot be isometrically embedded in a fiat 3D 
space. 

Eor odd multiple moment ci, the shadow is illustrated 
in figures [7][^ Here we observe an “eyebrow” structure 
for large values of multiple moments (e.g. ci = 1/800), 
reminiscent of those seen for the shadow of two merging 
black holes [35]. Eor an observer in the equatorial plane 
the shadow looks like an incomplete ellipse that is flat¬ 
tened on one side; furthermore the shape of the shadow is 
mirror-reflected with respect to the axis upon chang¬ 
ing the sign of the multipole moment from positive to 
negative. Eor an observer in the Oq = tt/A plane, we see 
the negative odd multiple moment to have a shape sim¬ 
ilar to that of the positive even multiple moment in the 
same plane; however the width of the shapes are smaller 
than the ones for even multiple moment. 
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VI. SUMMARY 

In many astrophysical situations, such as a black hole 
in a binary system, the black hole is distorted. Light rays 
emitted from sources behind the black hole will either 
be absorbed by it or escape to an observer located at 
some finite distance {vo^Oq) from the black hole. If a 
photon gets emitted from a point located on a sphere of 
radius re, the boundary of the set of rays reaching this 
observer bounds the ‘local shadow’ of the black hole for 
this observer. 

We have introduced this notion of local shadow in pre¬ 
cise terms and computed it for a distorted Schwarzschild 
black hole by tracing back the trajectory of these pho¬ 
tons. A null ray traced back from the point (ro, 6>o), with 
initial parameters will be part of the local shadow 

if this null ray is absorbed by the black hole; if it reaches 
a radius rg (after propagating in the space-time) it is not 
part of “local shadow” of the black hole. The parame¬ 
ters are related to the initial velocity integral 

of motion and the position of the observer. 

For a quadrupole distortion, C 2 7 ^ 0 we found that 
the shadow was deformed from a circle to an ellipse. As 
the magnitude of the quadrupole moment C 2 increases, 
the shadow gets increasingly deformed . Rather unex¬ 


pectedly, we found that positive values of C 2 the shadow 
is prolate whereas the horizon is oblate; for C 2 < 0 
the reverse occurs, with a prolate horizon having an 
oblate shadow. We have shown analytically why this phe¬ 
nomenon occurs for a local observer. However we expect 
that an observer at infinity will see the apparent shape of 
an oblate black hole (i.e., its shadow) to be either oblate 
or a circle of greater radius than for the undistorted cas^ 
The horizon is very rigid, i.e., it is harder to deform the 
horizon and easier to have a deformed shadow. In this 
paper, we don’t consider large enough distortions which 
distort the horizon effectively, since we do not want the 
distorting potentials to dominate over the potential of 
the black hole. 


ACKNOWLEDGMENTS 

The authors are grateful to the Natural Sciences and 
Engineering Research Council of Canada for financial 
support. The authors S. A. and C. T. also acknowl¬ 
edge support by the DEG Research Training Group 1620 
“Models of Gravity”. The authors would like to acknowl¬ 
edge Don N. Page for a useful comment. 


[1] S. S. Doeleman et al., Event-horizon-scale structure in the 
supermassive black hole candidate at the Galactic Centre^ 
Nature 455, (2008) 78 [astro-ph/0809.2442] 

[2] J. L. Synge, The escape of photons from gravitationally 
intense stars, Mon. Not. R. Astron. Soc. 131 (1966) 463. 

[3] J. M. Bardeen, Timelike and null geodesics in the Kerr 
metric, in: Black Holes-les astres occlus, edited by C. De- 
Witt and B. S. DeWitt (Gordon and Breach, New York, 
1973) p. 215 

[4] Z. Li and C. Bambi, Measuring the Kerr spin parameter 
of regular black holes from their shadow, J. Cosmol. As- 
tropart. Phys. JCAP 1401:041 (2014) arXiv:I309.I606 

[5] K. Hioki and K.-I. Maeda, Measurement of the Kerr Spin 
Parameter by Observation of a Compact Object’s Shadow, 
Phys. Rev. D 80 (2009) 024042 arXiv:0904.3575 

[6] T. Johannsen and D. Psaltis, Testing the No-Hair Theo¬ 
rem with Observations in the Electromagnetic Spectrum. 
I. Black-Hole Images, ApJ. 716 (2010) 187. 

[7] A. de Vries, The apparent shape of a rotating charged 
black hole, closed photon orbits and the bifurcation set 
A 4 , Class. Quantum Grav. 17 (2000) 123 

[8] P. Bakala, P. Cermak, S. Hledik, et al., Extreme gravita¬ 
tional lensing in vicinity of Schwarzs child-de Sitter black 
holes, Gen. Eur. J. Phys. 5 (2007) 599 

[9] C. Bambi and N. Yoshida, Shape and position of the 
shadow in the 6 = 2 Tomimatsu-Sato space-time. Class. 
Quant. Grav. 27 (2010) 205006 arXiv: 1004.3149 


^ This oblate/circle ambiguity is due to the fact that we cauuot 
deduce the value of at iufiuity aualytically. 


[10] L. Amarilla, E. F. Eiroa, and G. Giribet, Null geodesics 
and shadow of a rotating black hole in extended Chern- 
Simons modified gravity, Phys. Rev. D 81 (2010) 124045 
[gr-qc/I005.0607] 

[11] J. Schee and Z. Stuchlik, Optical phenomena in the field 
of braneworld Kerr black holes, Int. J. Mod. Phys. D 18 
(2009) 983 

[12] L. Amarilla and E. F. Eiroa, Shadow of a rotating 
braneworld black hole, Phys. Rev. D 85 (2012) 064019 
[gr-qc/III2.6349] 

[13] L. Amarilla and E. F. Eiroa, Shadow of a Kaluza-Klein 
rotating dilaton black hole, Phys. Rev. D 87 (2013) 
044057 arXiv:I30I.0532 

[14] A. Abdujabbarov, F. Atamurotov, Y. Kucukakca, B. 
Ahmedov, and U. Camci, Shadow of Kerr-Taub-NUT 
black hole, Astrophys. Space Sci. 344 (2012) 429 
arXiv:I2I2.4949 

[15] A. Grenzebach, V. Perlick, C. Lammerzahl, Photon re¬ 
gions and shadows of Kerr-Newman-NUT black holes 
with a cosmological constant, Phys. Rev. D 89 (2014) 
124004 arXiv: 1403.5234 

[16] A. Yumoto, D. Nitta, T. Chiba, and N. Sugiyama, Shad¬ 
ows of multi-black holes: Analytic exploration, Phys. Rev. 
D 86 (2012) I0300I arXiv: 1208.0635 

[17] J.W. Moffat, Black holes in modified gravity (MOG), 
arXiv:I4I2.5424 

[18] U. Papnoi, F. Atamurotov, S. G. Ghosh, B. Ahmedov, 
Shadow of five-dimensional rotating Myers-Perry black 
hole, Phys. Rev. D 90 (2014) 2, 024073 arXiv: 1407.0834 

[19] H. Falcke, F. Melia, and E. Agol, Viewing the shadow of 
the black Hole at the Galactic Genter, ApJ 528 (2000) 
LI3 





12 


[20] C. T. Cunningham and J. M. Bardeen, The optical ap¬ 
pearance of a star orbiting an extreme Kerr black hole^ 
ApJ 183 (1973) 237 

[21] J.-P. Luminet, Image of a spherical black hole with thin 
accretion disk, Astron. Astrophys. 75 (1979) 228 

[22] J. Dexter, E. Agol, P. C. Fragile, and J. C. McKinney, 
Radiative models of Sagittarius A * and M87 from rela¬ 
tivistic MHD simulations, J. Phys.:Con. Ser. 372 (2012) 
012023 [arXiv: 1202.0348] 

[23] M. Moscibrodzka, H. Shiokawa, C. F. Gammie, and J. C. 
Dolence, The Galactic Center weather forcast, ApJ Lett. 
752 (2012) LI arXiv:1204.1371 

[24] J. Dexter and P. C. Fragile, Tilted black hole accretion 
disc models of Sagittarius A time-variable millimetre 
to near-infrared emission, Mon. Not. R. Astron. Soc. 432 
arXiv: 1204.4454 (2013) 2252 

[25] R. Geroch and J. B. Hartle, Distorted black holes, J. 
Math. Phys. 23 (1982) 680 

[26] W. Israel and K. A. Khan, Collinear particles and Bondi 
dipoles in general relativity, Nuovo Cimento 33 (1964) 
331. 

[27] S. Chandrasekhar, The Mathematical Theory of Black 
Holes, Clarendon Press, Oxford University Press, New 
York, Oxford, (1983), pg. 583 

[28] A. Tomimatsu, Distorted rotating black holes, Phys. Lett. 


A 103 (1984) 374. 

[29] N. Breton, T. Denisova, and V. Manko, A Kerr black 
hole in the external gravitational field, Phys. Lett. A 230 
(1997) 7. 

[30] P. Peters, Toroidal black holes?, J. Math. Phys. 20 (1979) 
1481. 

[31] B. Xanthopoulos, Local toroidal black holes that are static 
and axisymmetric, Proc. R. Soc. Lond. A 388 (1983) 117. 

[32] S. Fairhurst and B. Krishnan, Distorted black holes 
with charge, Int. J. Mod. Phys. D 10 (2001) 691 gr- 
qc/0010088 

[33] S. Abdolrahimi, V. P. Frolov, A. A. Shoom, Interior of 
a Charged Distorted Black Hole, Phys. Rev. D 80 (2009) 
024011. 

[34] S. Abdolrahimi, A. A. Shoom, D. N. Page, Distorted 5- 
dimensional vacuum black hole, Phys. Rev. D 82 (2010) 
124039. 

[35] R. Cowen, Black-hole mergers cast kaleidoscope of shad¬ 
ows, Nature (2014) 16283. 

[36] M. Berger, A Panoramic View of Riemannian Geometry, 
Springer-Verlag, Berlin-Heidelberg-New York, (2003). 

[37] E. F. Taylor and J. A. Wheeler, Exploring black holes, 
Introduction to general relativity, Addison Wesley Long¬ 
man, Inc. USA, 2000, pg. 5-11. 


